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ABSTRACT 

The Dark Energy Camera has captured a large set of images as part of Science Verifi¬ 
cation (SV) for the Dark Energy Survey. The SV footprint covers a large portion of the 
outer Large Magellanic Cloud (LMC), providing photometry 1.5 magnitudes fainter 
than the main sequence turn-off of the oldest LMC stellar population. We derive geo¬ 
metrical and structural parameters for various stellar populations in the LMC disk. Eor 
the distribution of all LMC stars, we find an inclination of i = —38.14° ± 0.08° (near 
side in the North) and a position angle for the line of nodes of Oq = 129.51° ±0.17°. We 
find that stars younger than ^ 4 Gyr are more centrally concentrated than older stars. 
Eitting a projected exponential disk shows that the scale radius of the old populations 
is R>4Gyr = 1-41 ± 0.01 kpc, while the younger population has R^AGyr = 0.72 ± 0.01 
kpc. However, the spatial distribution of the younger population deviates significantly 
from the projected exponential disk model. The distribution of old stars suggests a 
large truncation radius ofRt = 13.5±0.8 kpc. If this truncation is dominated by the 
tidal field of the Galaxy, we find that the LMG is 24lg times less massive than the 
encircled Galactic mass. By measuring the Red Glump peak magnitude and comparing 
with the best-fit LMG disk model, we find that the LMG disk is warped and thicker 
in the outer regions north of the LMG centre. Our findings may either be interpreted 
as a warped and flared disk in the LMG outskirts, or as evidence of a spheroidal halo 
component. 

Key words: galaxies: Magellanic Glouds; galaxies: stellar content; stars: statistics 


1 INTRODUCTION 

The Milky Way (MW) satellite system offers a variety of ex¬ 
amples of dwarf galaxies. Most of its members are essentially 
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gas-free and contain mainly old stars (McConnachie 2012). 
The evolution of these systems is closely related to the for¬ 
mation of the Galaxy and the process of mass assembly of 
the large scale structures in the Universe (Klypin et al. 1999; 
Moore et al. 1999; Stewart et al. 2008). On the other hand, 
the Large and Small Magellanic Clouds (LMC and SMC, 
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respectively) are the closest low-mass, gas-rich (Grcevich & 
Putman 2009) interacting systems. The main features trac¬ 
ing the interaction history of the clouds are the HI Mag¬ 
ellanic Stream (Mathewson, Cleary & Murray 1974), and 
Bridge (Hindman, Kerr & McGee 1963). A counterpart of 
the Stream was found and named the Leading Arm (Putman 
et al. 1998). The formation of these structures is a subject 
of great debate. Recent simulations favour a scenario where 
the Stream, Bridge, and Leading Arm are remains from the 
close interaction between the LMC and SMC before falling 
into the MW potential (Besla et al. 2012; Kallivayalil et al. 
2013). 

The Clouds’ star formation history (SFH) also reflects 
their close interaction history. It is possible to identify multi¬ 
ple periods of enhanced star formation that are arguably cor¬ 
related with close encounter between the Clouds (Meschin 
et al. 2014; Rubele et al. 2012; Javiel, Santiago & Kerber 
2005; Holtzman et al. 1999). Evidence of such events are im¬ 
printed in the stellar population which, due to the proximity 
of the Magellanic System, is resolved into individual stars 
with medium sized ground-based telescopes. Enhanced star 
formation is also demonstrated by the extensive star clus¬ 
ter system throughout the Magellanic Clouds (MC). The 
clusters in this system span a very broad range in age and 
metallicity (Kerber & Santiago 2009). Evidence of an age- 
gap (Jensen, Mould & Reid 1988) may support a relationship 
between the formation/disruption rate of star clusters and 
the intergalactic interaction history. In this sense the star 
clusters in the LMC may give hints to how the intergalac¬ 
tic interaction affects the evolution of star cluster systems 
(Renaud & Gieles 2013). 

Despite being the nearest interacting system of galax¬ 
ies, the Magellanic System still has only a small angular 
fraction observed to the photometric depth of its old main 
sequence turn-off (MSTO). Deep observations suggest that 
the LMC stellar populations extend beyond an angular dis¬ 
tance of 15° from its centre (Majewski et al. 2009). There 
is also kinematic evidence for a dynamically warm stellar 
component consistent with a halo (Minniti et al. 2003) that 
has its major axis oriented with the disk (Alves 2004). 

Very few studies are available in the outskirts of the 
LMC. Weinberg V Nikolaev (2001) report an exponential 
scale length of Rs ~ 1.4 kpc with no significant distinction 
between a young and old disk; however, their sample is from 
the 2 Micron All-Sky Survey (2MASS; Skrutskie et al. 2006), 
which is very shallow and does not allow for a clear age 
selection. Saha et al. (2010) report a smaller scale length of 
Rs = 1.15 kpc based on an optical survey. They argue in 
favour of a truncation radius of Rt ~ 14 kpc. Their sample 
is limited to only a few fields and their analysis does not use 
age selected stellar samples. 

A new generation of photometric surveys is now com¬ 
ing to the Southern Hemisphere, allowing for the first time 
a complete view of the Magellanic System. One of them 
is the Dark Energy Survey (DES), which will observe the 
outskirts of the LMC, SMC, and most of the Magellanic 
Stream over the course of 5 years. DES is a photometric 
survey with the primary goal of measuring the dark energy 
equation of state. To achieve this goal the survey will em¬ 
ploy four independent cosmological probes: galaxy clusters, 
baryon acoustic oscillations, weak leasing, and type la su¬ 
pernovae (Elaugher 2005). The total survey area is ~ 5000 


deg^ reaching a magnitude limit of z ~ 24. The photomet¬ 
ric system adopted for DES comprises the filters r, z, z, 
which are similar to the ones used in the Sloan Digital Sky 
Survey (SDSS; Eukugita et al. 1996), with the addition of 
the Y passband, which provides synergy with the VISTA 
Hemisphere Survey (McMahon 2012; Cioni et al. 2011). The 
DES footprint will overlap with several other surveys in the 
southern hemisphere allowing a multi wavelength approach 
to various astrophysical problems. 

An early release of DES Science Verification (SV) data 
was made available to the DES collaboration recently. The 
data cover ~200 deg^ of the southern sky and sample regions 
as close as 4° North from the LMC centre. The photometric 
catalogue from this release reaches ~ 1.5 mag fainter than 
the old MSTO of the LMC (which is at ^ ~ 22), allowing 
for a detailed study of the resolved stellar population of this 
galaxy. 

In this paper we use the DES SV data to study the LMC 
geometry and density profile as traced by stellar components 
with a characteristic age range. We model the distribution of 
stars using a simple projected exponential disk and perform 
a formal fit. We discuss the presence of a possible trunca¬ 
tion radius and its implication for the LMC mass. As an 
alternative probe of the LMC geometry we use Red Clump 
(RC) stars as a distance indicator and a ruler for the LMC 
thickness. This paper is organized as follows. In Section 2 
we present a brief introduction to the DES SV data. In Sec¬ 
tion 3 we discuss the quality of the photometry and address 
issues due to completeness and survey coverage. In Section 
4 we describe the disk model and the fitting procedure used 
to find the geometrical parameters of the LMC. Section 5 
shows our efforts to use the RC as a distance and thick¬ 
ness estimator. In Section 6 we summarize and discuss the 
implications of the results found in this paper. 


2 DECAM AND DES SCIENCE 

VERIFICATION 

The Dark Energy Camera (DECam) (Elaugher et al. 2010) 
was constructed in order to carry out the Dark Energy 
Survey. This instrument has a focal plane comprised of 74 
CCDs: 62 2kx4k CCDs dedicated to science imaging and 12 
2k x2k CCDs for guiding, focus, and alignment. The camera 
is also equipped with a five element optical corrector and a 
sophisticated cryogenic cooling system. DECam is installed 
at the prime focus of the Cerro Tololo Inter-American Ob¬ 
servatory (CTIO) 4 meter Blanco telescope. In this config¬ 
uration, DECam has a 2.2 degree-wide field-of-view (EoV) 
and a central pixel scale of 0.263 typical site con¬ 

ditions, the Blanco telescope plus DECam yield an imaging 
point spread function (PSE) with full-width half maximum 
(EWHM) of 0.9'^ which is adequately sampled by the pixel 
scale. 

DECam was commissioned in September 2012 and be¬ 
gan operations in November 2012, with the DES SV cam¬ 
paign covering a period of three months. The DES SV data 
is intended to test capabilities of the camera, the data trans¬ 
fer infrastructure, and the data processing. All images taken 
during SV are public; however, the catalogues generated by 
the collaboration are proprietary. 
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2.1 Data reduction 

The DES data management (DESDM) team was responsible 
for the reduction of the SV images. Here we will give a brief 
description of the reduction process. Eor a complete descrip¬ 
tion we refer to Sevilla et al. (2011), Desai et al. (2012), and 
Mohr et al. (2012). 

The DESDM data reduction pipeline consists of the fol¬ 
lowing steps: 

Image detrending: this step includes correction for 
crosstalk between CCD amplifier electronics, bias level cor¬ 
rection, correction for pixel-to-pixel sensitivity variation 
(flat-fielding), as well as corrections for non-linearity, fring¬ 
ing, pupil, and illumination. 

Astrometric calibration: in this step, bright known 
stars are identified in a given image using S Extractor 
(Bertin & Arnouts 1996). The position of these stars in the 
focal plane is used to find the astrometric solution with the 
aid of the software SCAMP (Bertin 2006) through compar¬ 
ison to UCAC-4 (Zacharias et al. 2013). 

Nightly photometric calibration: several reference 
stars are observed each night, and a photometric equation 
is derived for that night. This equation takes into account 
a zero point, a colour term, and an airmass term for each 
of the DECam science CCDs. 

Global calibration: relative photometric calibration of 
the DES survey area is done with repeated observations 
of the same star in overlapping DECam exposures using a 
method similar to that described in Glazebrook et al. (1994). 
The relative magnitudes are anchored to a small set of abso¬ 
lutely calibrated reference stars within the same area called 
tertiary standards which come from observations on photo¬ 
metric nights. They are then calibrated relative to known 
equatorial belt standards observed on the same night with 
the same filter set (Tucker et al. 2007). The current accuracy 
of the relative and absolute systems is a few percent, and is 
expected to improve as the DES survey covers larger con¬ 
tiguous areas. DES calibration was checked using the Star 
Locus Regression technique (Kelly et al. 2014). This anal¬ 
ysis revealed that zero-point offsets are typically less than 
0.05 mag in r, z, and z across the DES footprint observed 
so far. 

Coaddition: in order to increase the signal-to-noise ra¬ 
tio, exposures are combined. This has the advantage of mit¬ 
igating transient objects, such as cosmic rays and satellite 
trails. This step requires the placement of the images in a 
common reference projection. The software SWarp (Bertin 
et al. 2002) was used and the coaddition process was done 
in segments of the sky called tiles. At this stage the flux is 
corrected according to the photometric calibration described 
above. 

Cataloging: for each coadded tile source detection and 
model-fitting photometry is performed using PSEEx and 
SExtractor (Bertin 2011; Bertin & Arnouts 1996) on a 
combined r, z, and z image. Object fluxes and many other 
characteristics are calculated in the individual g,r,i,z,Y 
frames. These catalogues are ingested into a high perfor¬ 
mance database system. 

The final catalogue is available for the DES collabo¬ 
ration through a database query client. There are approxi¬ 
mately 900 parameters measured for each source identified 


Table 1. Here we show a subset of the SExtractor parameters 
measured for the coadded DECam SV data. We also show any 
selection criteria that were made. 


Parameter name Selection 

RA 
DEC 

MAG_AUT0_* 

MAGERR_AUT0_* 

SPREAD_M0DEL_* 

FLAGS.* 

The * symbol refers to all passbands available. 


by SExtractor. In Table 1 we list a few parameters relevant 
to this work. In the same table we also list any comments 
about the parameter and the quality cuts applied. 


3 THE SV DATA 

By the end of the SV campaign, DECam had obtained im¬ 
ages in five passbands of a region with roughly 300 deg^, 
for which 200 deg^ are contiguous. The contiguous region 
covers the northern outskirts of the LMC. This region over¬ 
laps an eastern portion of the South Pole Telescope (SPT) 
footprint (Carlstrom et al. 2011). Hence we call it SPT-E 
for simplicity. 

In this section we will discuss several aspects of the SV 
data and how it is suitable for the analysis we propose. 


3.1 Photometry 

Currently, the deepest and most homogeneous magni¬ 
tude measurement resulting from the DESDM pipeline is 
MAG_AUT0, which is computed using the flux inside a Kron 
radius (Kron 1980). The Kron radius is dependent on how 
extended a source is, hence the aperture is variable, but it is 
essentially the same for all stars, given that there are no sig¬ 
nificant spatial dependence in the image quality, which is the 
case for DES observations. Thus, MAG_AUT0 is roughly equiv¬ 
alent to a simple aperture magnitude but less sensitive to 
seeing variations. We adopt the following notation through¬ 
out the paper: g is the MAG_AUT0 magnitude measured for the 
g passband. This also applies for the magnitudes measured 
in r, z, z, Y. 

The raw DESDM catalogue has approximately 10® 
sources with g ^ 24.6. This list includes spurious detections 
like satellite trails, star wings, cosmic rays, etc. To exclude 
such detections from our catalogue we adopt a simple cut in 
the FLAGS parameter. We select only sources that simultane¬ 
ously have FLAGS_G and FLAGS_R ^ 3, which selects objects 
that are not saturated and do not contain any bad pixel. The 
FLAGS code is the same as the one adopted by SExtractor. 
The number of sources left after this cut is ~ 9 x 10^. 

We check the stability of the photometric calibration 
provided by DESDM by comparing the RC peak colour at 
different points of the SPT-E that contain LMC stars. Eor 
all bands we found a maximum scatter of 0.02 mags around 
the RC peak colour. 


^ 26.0 in g and r 

|SPREAD_M0DEL_I| < 0.002 

^ 3 in g and r 
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Table 2. Some general photometric system information: the cen¬ 
tral wavelength (Ac) in nanometres, the extinction as a fraction 
of the extinction in the Johnson V passband (assuming Ry = 3.1 
and a Cardelli, Clayton &: Mathis (1989) extinction curve for the 
MW), and the value of the magnitude corresponding to a typical 
signal-to-noise ratio of 10. In the last three columns the coeffi¬ 
cients of the best-fit error model are shown. 


Filter 

Ac 

(nm) 

Ax/Av 

^10% 

(a, 6, c) 

g 

479 

1.199 

23.94 

(0.001, 26.41, 1.25) 

r 

641 

0.837 

23.76 

(0.001, 26.34, 1.27) 

i 

781 

0.635 

22.75 

(0.003, 25.52, 1.34) 

z 

924 

0.462 

22.03 

(0.003, 24.75, 1.43) 

Y 

1008 

0.400 

20.50 

(0.009, 23.45, 1.40) 


The following model was adopted to describe the pho¬ 
tometric uncertainties from the DECam SV data. 

a[mag) = a + exp I - 1 (1) 

By fitting the above equation to a sample of 0.01% ran¬ 
domly chosen stars from the SPT-E region we find the error 
curves shown in Eigure 1. The coefficients for each curve are 
given on Table 2, where we also show the 10% uncertainty 
magnitude for each band. Along with these values we also 
give some basic information about the filter central wave¬ 
length and extinction coefficients at those wavelengths from 
Cardelli, Clayton & Mathis (1989). 

3.2 Star/Galaxy Separation 

The DES collaboration has generated several datasets for 
validating and testing the DESDM system (Mohr et al. 
2012). These datasets were simulations of the actual obser¬ 
vations. To create these mock observations, input catalogues 
of artificial galaxies and stars were generated (Rossetto et al. 
2011; Balbinot et al. 2012; Busha 2014). The mock obser¬ 
vations included a time varying seeing, realistic shapes for 
the galaxies, and variations to the focal plane of DECam. 
These images were fed into the DESDM reduction pipeline 
and an output catalogue was generated, which was then re¬ 
leased to the collaboration to perform tests of their scientific 
algorithms. 

Tests with a set of such simulations, called Data Chal¬ 
lenge 5 (DC5) by the collaboration, have been carried 
out in order to assess how different star-galaxy classifica¬ 
tion parameters perform. A summary of this comparison 
is shown in Rossetto et al. (2011). The main conclusion is 
that SPREAD_M0DEL (Desai et al. 2012; Bony et al. 2013) 
performs better in terms of purity and completeness than 
other typically employed classifiers such as CLASS_STAR and 
FLUX_RADIUS (Bertin & Arnouts 1996). Later, Soumagnac 
et al. (2013) developed a more sophisticated star-galaxy sep¬ 
aration algorithm; however, this algorithm must be trained 
on a data set where true stars and galaxies are known. Im¬ 
plementation of such methods are being considered collabo¬ 
ration wide. Eor simplicity, we choose to use SPREAD_M0DEL 
measured in the i band as the star-galaxy separator. The 
cut-off criterion for selecting stars is |SPREAD_M0DEL_l| ^ 
0.002. Erom test on DC5, this cut is found to correspond to 
a simultaneous stellar completeness and purity of ~ 80% for 






Figure 1. Magnitude versus uncertainty for a sample of 0.01% 
of the total number of stars in the DECam SV region. The solid 
curve shows the best-fit error model. The coefficients are given in 
Table 2. The dashed line shows the 10% uncertainty magnitude. 

objects with g < 23 (Rossetto et al. 2011). Erom this point 
further we call objects that match this cut-off criteria stars. 
It is worth mentioning that at ^ = 23 we expect ~ 20% 
of contamination from galaxies in our sample. However, the 
large scale distribution of these objects is homogeneous and 
is unlikely to significantly affect the findings of this paper. 

3.3 Completeness 

To independently assess the completeness of the DESDM 
catalog we conducted a few experiments using DAOPHOT 
(Stetson 1987). This photometry code is known to perform 
very well in extremely crowded regions such as the cores of 
globular clusters (Balbinot et al. 2009). Hence, at the typical 
density of LMC field stars it should yield a fairly complete 
catalogue. This DAOPHOT catalogue may be compared 
to the one produced by DESDM as an approximation to a 
complete catalogue and giving an estimate of the complete¬ 
ness as function of magnitude. This is not the most accurate 
approach to this problem; however, it is much simpler than 
performing artificial star experiments across several hundred 
squared degrees. 

We performed DAOPHOT photometry on 51 fields, 50 
of which contained a LMC star cluster in its centre. The 
51®^ field was selected as far away as available in the SPT-E 
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data in order to sample a region where there are few or no 
LMC stars. The fields selected are subregions of the coadded 
images encompassing 6.75' x 6.75' each, with the exception 
of the 51®^ field, which has 18' x 36'. The first 50 fields 
were selected in order to assess the completeness not only 
in regions with a density of stars typical of the LMC, but 
also with varying density, such as the inner regions of a 
star cluster. A broader discussion about this subject will be 
presented in a future paper. The 51®^ field was selected in 
order to compare the performance of the DESDM reduction 
pipeline in a region with little or no crowding. The position 
of each of the all 51 fields is marked on Figure 2. On average 
DESDM and DAOPHOT photometry agree within 0.02 in 
g and r. 

To compare the number of stars as a function of magni¬ 
tude, we cut both the DESDM and DAOPHOT catalogues 
at the 3% photometric uncertainty level. This cut happens 
at g ~ 23.5. To separate stars from galaxies in DAOPHOT 
we adopted a cut in the sharpness parameter which behaves 
similarly to SPREAD_M0DEL. 

In Figure 3 top panel we show the ratio between the 
number of stars detected by DESDM and DAOPHOT 
{Ndes/Ndaophot) averaged in magnitude bins for the 51 
fields discussed above. The solid line shows the average value 
at a given magnitude bin. The shaded region shows the stan¬ 
dard deviation. In the two bottom panels in the same figure 
we show the value of Ndes /Ndaophot as a function of 
angular separation to the LMC centre. We show this for 
two magnitude bins: 21.8 < ^ 22.6 and 22.6 < g ^ 23.3 

(also indicated in the figure). From Figure 3 we notice that 
the average completeness drops as a function of increasing 
magnitude, as expected. The completeness remains roughly 
constant as a function of angular separation to the LMC 
centre with a spread of the same order as the one shown 
in the shaded region of the top panel. We conclude that 
the completeness has little spatial dependence at least up to 
magnitudes as faint as ^ ~ 23. 

Despite the good performance in moderately crowded 
fields (i.e. LMC field population), we observed that in the 
inner regions of star clusters the DESDM sample is highly 
incomplete reaching less than 40% completeness for ^ = 19 
and dropping to zero for fainter magnitudes. This will be 
further explored in future publications. 

We note that where the completeness is most important 
to our results regarding truncation radii and scale lengths, 
i.e. at distances farther away from the center of the LMC, 
the star density is low enough that crowding effects do not 
affect the sample selection as shown in Figure 3. 

From the completeness analysis described above we 
conclude that in the DESDM catalogue completeness does 
not depend strongly on source density, except in extremely 
crowded environments such as the central parts of stars clus¬ 
ters. From visual inspection and through comparison with 
state-of-the-art crowded field photometry methods we find 
no strong evidence for large-scale variations in the stellar 
completeness in the SPT-E region, thus allowing us to prop¬ 
erly access the spatial distribution of stars across its foot¬ 
print, either from the LMC or from the MW. 


3.4 Mangle masks 

We use the Mangle software (Swanson et al. 2008) to track 
the survey coverage and limiting magnitude. Along with the 
release of the DESDM catalogues, a set of Mangle masks 
was also provided. These masks contain information about 
the limiting magnitude in each patch of the sky, as well as 
other key information about data quality. The magnitude 
limit is computed as the detection limit of a point source 
with signal-to-noise ratio S/N = 10. The noise is estimated 
from the variance and level of the sky on a patch of the sky. 
To measure the flux an aperture of l" radius is used; this 
yields a magnitude measurement called MAG_APER_4. For a 
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Figure 2. Top panel: Gnomonic projection of the number of 
2MASS stars that are consistent with AGB and RGB stars ac¬ 
cording to Yang et al. (2007). The projection is centred in the 
LMG centre. Part of the SMG is visible in the lower right cor¬ 
ner. The solid contour shows an approximate footprint of the 
southern part of the SPT-E region of the DES SV data. Bottom 
panel: Gnomonic projection of a Nsi ( i ^= 4:096 HEALpix map of 
the g magnitude limit Mangle mask for the SPT-E region quoted 
above. Holes and unobserved regions are masked and shown in 
white. Red circles mark the position of the fields selected for as¬ 
sessing the survey completeness. 
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Figure 3. Top panel: we show the ratio of stars detected by 
DESDM and DAOPHOT as a function of g magnitude. The solid 
line shows the average for the 51 fields analysed. The shaded 
region shows the standard deviation for the 51 fields. The dashed 
vertical line shows the faintest magnitude limit used in this work. 
Bottom panels: black circles show the same ratio as in the top 
panel but now as a function of angular separation to the LMC 
centre for two magnitude bins (indicated in the captions). The 
error bars are computed using the Poisson uncertainty. 

complete description of the masks we refer to Swanson et al. 

( 2012 ). 

In Figure 2 we show the number of AGB and RGB 
stars found in 2MASS according to Yang et al. (2007) with 
the southern limits of the SPT-E region overplotted (top 
panel). We also show an approximation of the magnitude 
limit Mangle mask for the SPT-E region (bottom panel). It 
is an approximation in the sense that it uses the value of the 
Mangle mask at the central position of each pixel in the sky. 
The gray scale represents the magnitude limit at each point 
of the footprint. Holes caused by bright stars or other im¬ 
age imperfections are displayed as well. These regions have 
a magnitude limit of zero, hence going out of the gray scale 
range. The white color in the figure represents regions that 
were not observed. The figure uses a HEALPIX pixelization 
(Gorski et al. 2005) with Nside=^096 and a Gnomonic pro¬ 
jection centered at (a,(5) = (75°,—55°). 

The approximate mask discussed above allows us to deal 
with the full SPT-E mask in a much faster way and at the 
same time to provide a very good approximation of the gen¬ 
eral properties of the survey, such as coverage and magnitude 
limits. The coverage mask is simply a HEALPIX map with 
the value of 1 for pixels that were observed and 0 for those 
that were not or were masked for some reason. The coverage 
mask originally did not contain holes due to star clusters. 
To mask these regions we conducted a visual search for star 


Figure 4. The Hess diagrams in each set of colour and magnitude 
available in the DES photometric system. In the top left panel 
we also show the colour-magnitude polygons used to select stars 
consistent with a young/intermediate (dot dashed polygon) and 
old (solid polygon) stellar population. The boxes (dashed) used 
for the selection of RC stars are shown in all panels. 


clusters in the SPT-E and used the position of each cluster 
to add a hole in the mask (i.e. a region with coverage equals 
to zero). 

The approximated masks were used to remove stars re¬ 
siding in regions where the following criteria were met: (i) 
a limiting g or r magnitude brighter than 23; (ii) coverage 
value equal to 0. The additional trimming for zero mask val¬ 
ues was necessary because the original catalogue was limited 
by the exact masks, not the approximated ones. 

The above section outlines the process of trimming the 
catalogue to select objects that are likely stars and to keep 
only regions with deep photometry that were not affected 
by artifacts in the survey. This catalogue, as well as the 
mask associated with it, will be used in the remainder of 
this paper. 


4 THE LMC GEOMETRY 

Glassically, the LMG is classified as an Irregular Dwarf 
Galaxy, although it has several major components of a spiral 
galaxy such as a disk and a bar. It is, perhaps, more appro¬ 
priate to classify this galaxy as a highly perturbed spiral 
galaxy. Its morphology departs so much from a classical Ir¬ 
regular Dwarf that it has been established that the LMG is 
the prototype of a class of dwarf galaxies called Magellanic 
Irregular (de Vaucouleurs & Ereeman 1972). These galaxies 
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are characterized by being gas-rich, one-armed spirals with 
off-centre bars. 

The SMC is the closest neighbour to the LMC. Together 
they form the Magellanic System. Recent dynamical mod¬ 
elling (Besla et al. 2012) and high precision tangential veloc¬ 
ity measurements (Kallivayalil et al. 2013) point to the need 
for updated thinking with regard to the origins of the Mag¬ 
ellanic System. The centre of mass spatial velocity of these 
galaxies is very close to the escape velocity of the MW, thus 
suggesting that the system is not gravitationally bound to 
the MW. The same models also predict the formation of the 
Magellanic Bridge and Stream as a result of the interaction 
of the LMC with the SMC, generating long arms of debris 
in the same fashion as the Antennae system. 

Despite the growing number of simulations and high 
precision velocity measurements, a few aspects of the LMC 
geometry, such as the presence of a spheroidal halo (Ma- 
jewski et al. 2009), and the warping and flaring of the disk 
(Subramaniam & Subramanian 2009) remain uncertain to 
some degree. New large area surveys in the Southern Hemi¬ 
sphere are beginning to shed light on these uncertainties. 

Here we study the LMC disk geometry using a very 
simple approach. We try to model its stellar density using a 
circular exponential disk. This disk is inclined relative to the 
sky plane by the angle i. To compute the expected number 
of stars p as a function of a and 5 we use the transformations 
found in Weinberg & Nikolaev (2001). For the sake of clarity, 
we give the expression for the heliocentric distance t to a 
given point of the disk with coordinates (a, (5). 

t((a, ^) = —Rlmc cos i x {cos S sin (a — ao) sin 0 sin i 

+ [sin S cos ^0 — cos S sin cos (a — ao)] cos 0 sin i 
— [cos S cos ^0 cos {a — ao) + sin S sin ^o] cos i}~^ (2) 

where (ao, ^o) is the central coordinate of the LMC, i is the 
disk inclination, Rlmc is the heliocentric distance to the 
LMC centre, and 0 the position angle (PA) of the minor 
axis. In the reference frame adopted here, the inclination i 
has negative values when the North side of the LMC is closer 
to us. In the reference frame adopted, the i angle is reversed 
to what is typically adopted in the literature. 

The density of stars is simply given by 

p(a, S) = pot{a, S)‘^exp{-R/Rs) + pBcfi^, S) (3) 

where R is the radial distance in the disk plane, Rs is the 
scale length of the exponential disk, pot(ao, ^o) is the central 
density, and pbg is the density of background/foreground 
stars. The function /(a, (5) is a third degree polynomial 
which takes into account the spatial variation of MW field 
stars. 

Five parameters are used to model the disk geometry: 
the LMC central coordinates ((ao,^o), its distance to the 
Sun (Rlmc), the inclination z, and the PA 0. Additionally, 
3 parameters describe the density of stars along this disk: 
the central density po, the scale radius Rg, and the density 
scale of background stars pbc- 

To simplify the problem and to better accommodate 
the fact that the observations are all on the northern side 
of the MCs, we make a few assumptions. The LMC cen¬ 
tre is kept fixed at the Nikolaev et al. (2004) value of 



Age [Gyr] 

Figure 5. Ratio of the number of young (old) to all stars as a 
function of age. The ratio was computed using a constant SFR 
and a Piatti & Geisler (2013) AMR. 


((ao,(^o) = (79.40°,—69.03°). The heliocentric distance to 
the LMC centre Rlmc = 49.9 kpc is adopted from the most 
recent review about the subject (de Grijs, Wicker & Bono 
2014). 

With the machinery to produce LMC disk models we 
proceed to a formal fit to the observed distribution of stars 
in the SPT-E region. This fit was performed on three sam¬ 
ples of stars. The selections were made based on PARSEC 
(Bressan et al. 2012) stellar evolution models. The first is 
made only of stars older than 4 Gyr, selected using a colour- 
magnitude cut shown in Eigure 4. The second sample se¬ 
lects only stars younger than 4 Gyr. The third sample is 
made of stars that fall within the limits 17 < p < 23 and 
—0.5 < g — r < I.O. We name these samples old, young, 
and all for simplicity. In Eigure 5 we show the ratio of old 
(young) and all stars for a simulated LMC stellar popula¬ 
tion that assumes a constant star formation rate (SER) and 
a Piatti & Geisler (2013) age-metallicity relation (AMR). 
Photometric errors were simulated according to Equation 
I and Table 2. We observe that the colour-magnitude cuts 
chosen are able to create two samples with virtually no age 
overlap. 

Using a TRILEGAL (Girardi et al. 2005) simulation, 
we modeled the expected MW stellar population in the SPT- 
E region. Assuming the same colour cuts as the ones used to 
select the three stellar populations quoted above, we fit the 
density of stars as a function of RA and Dec using a third 
order polynomial, f(a,S). The best-fit polynomial for each 
of the stellar populations cuts is later used to account for 
foreground contamination while fitting the LMC disk. This 
polynomial is scaled by pbg since the colour-magnitude cuts 
are not perfectly consistent with the theoretical prediction 
of TRILEGAL. 

To count stars in the sky plane we adopt the HEALPIX 
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pixelization scheme with Ngide — 512 which yields a con¬ 
stant pixel area across the sky of ~ 0.013 deg^. The maps 
built in this scheme for each stellar sample described above 
are shown in the left panel of Figures 6, 7, and 8. 

The best-fit disk model was obtained through a Markov 
Chain Monte Carlo (MCMC) technique. We used the code 
EMCEE (Foreman-Mackey et al. 2013) in its version 2.0.0. 
The test statistics chosen is a binned Poisson log-likelihood 
model (Dolphin 2002). We refer the reader to these authors 
for details on the MCMC and statistics used. This is the 
most appropriate test statistics for data-model comparison 
since we are dealing with counts that are subject to shot 
noise, especially in the outermost regions of the LMC. Our 
MCMC run uses a total of 30 walkers that make 1000 steps 
each for the hurn-in phase. After the hurn-in^ we let the 
walkers advance 5000 more steps each, sufficient to well sam¬ 
ple the parameter space and converge to the maxima. 

This fitting procedure was repeated for each stellar pop¬ 
ulation {all, old and young), yielding the parameters shown 
on Table 3. In this table we chose to present the PA of the 
line of nodes Oq = 0 90°. This is the quantity most often 

presented in the literature. 

The boxes chosen to select a given stellar population in 
the colour-magnitude diagram (CMD) assume that the stars 
are all at the same distance to us. This is not strictly the 
case for the LMC stars, which are spread across a disk plane 
inclined relative to the sky plane. To test how this distance 
spread affects our disk fit we correct the magnitude of each 
star using the best-fit disk models so as to bring the stars to 
a common distance. This common distance is chosen as the 
mean distance to the disk in the SPT-E. Using these distance 
corrected magnitudes we perform the CMD box selections 
again. For all the CMD based selections an increase of < 1% 
in the number of stars inside the CMD box is observed. Using 
this distance corrected sample we rerun the fitting procedure 
and find that the largest change in the parameters is a 5% 
increase in the inclination. The remaining disk parameters 
are more stable. We adopt the difference on the parameters 
found in this experiment as the systematic uncertainty in 
our parameter estimation. We show this variations as the 
uncertainties in parenthesis in Table 3. 

The left panel of Figure 6 a HEALPIX map in Gnomonic 
projection shows the logarithm of the number of stars for the 
all stellar population. Overploted we show the best-fit disk 
model as isodensity lines. On the right side of the figure we 
show the sampled probability distribution function for each 
parameter after marginalization. In light blue we mark the 
point of maximum probability density, which indicates the 
maximum likelihood solution. Notice that the samples are 
distributed in a fairly symmetric way around the maxima, 
yielding symmetric error bars. On Eigures 7 and 8 are similar 
to the figure described above, but showing different age- 
selected populations. 

We notice that the old population spreads out to de¬ 
clinations of ~ —55°, while the young population is much 
more abruptly truncated at ~ —60°. This points to different 
scale radii for these populations. In Table 3 we give a sum¬ 
mary of the parameters that best-fit each case. Here we see 
that the distributions have significantly different values of 
Rs while retaining similar values for the purely geometrical 
parameters Oq and i. 

The difference in the scale length is much more obvious 


in Eigure 9 where we show, for each stellar population, the 
average number of stars per HEALPIX pixel in bins of dis¬ 
tance along the LMC disk. We also show the best-fit disk 
model for each of the stellar populations, this model in¬ 
cludes the MW foreground stars (i.e. the term pBGf{a,5) 
in Equation 2). We notice that the old (triangles and dot- 
dashed line) profile is much more extended than the young 
one (crosses and dotted line). It is also remarkable that the 
LMC old density profile is well fit without the need for other 
components such as a spheroidal halo. However, the young 
profile is not very well described by the disk model. This 
could also account for the value found for its inclination, 
which is not in good agreement with the literature (van der 
Marel & Cioni 2001; Nikolaev et al. 2004; Rubele et al. 2012). 
In the same figure, the solid line shows the contribution of 
MW stars, we notice that the slope in the number of stars 
at the outskirts of the LMC can be accounted by for the 
spatial variation in the number of MW field stars. Due to 
the large number of stars in each radial bin, the Poisson un¬ 
certainty is very low. These uncertainties are only visible, 
as error bars, on the inset plot, where the outer parts of the 
old population profile is shown. 

We notice from Table 3 that the statistical uncertain¬ 
ties in the disk parameters are quite small. This is caused 
by the relatively poor description that our disk model gives 
of the LMC stellar population, especially the young popu¬ 
lation. In Eigure 9 this becomes apparent where we observe 
that the models deviate significantly from the observations, 
even when the uncertainties are considered. This can point 
to the case where there are other systematic effects that were 
not taken into account. One cause of such effects might be 
the spatial variation of the star-galaxy classifier. However, 
the classifier adopted here is very stable especially at bright 
magnitudes {g < 22.5), thus being improbable to be a sig¬ 
nificant source of error. Also, spatial changes in the com¬ 
pleteness of the survey are unlikely to cause a large change 
in the number of stars per HEALPIX pixel (see Section 3.3). 
Evidence indicates that the deviations from the fitted disk 
model are real features of the LMC structure. These features 
are more apparent for the young stellar population, which is 
most likely to hold signs of disk perturbations such as spiral 
arms. 

We use the old population profile to probe the total ex¬ 
tent of the LMC. We define the truncation radius as the ra¬ 
dius where the observed density profile reaches becomes in¬ 
distinguishable from the MW foreground population. In the 
inset plot of Eigure 9 we observe that for \og{R/kpc) >1.13 
the profile can be explained solely by the MW contribution. 
The uncertainty on this truncation radius (Rt) is taken as 
the size of the radial bin, which is 0.8 kpc. This yields a 
Rt = 13.5 ± 0.8 kpc. 

If we assume that the LMC luminous component is 
tidally truncated by the MW potential we can use the simple 
theoretical tidal radius formula (Binney & Tremaine 2008), 
given by: 

we find the following relation for the LMC mass (Mlmc) 
and the MW mass (Mmw) encircled within a radius equal 
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Population 

Age 

Gyr 

^0 

deg 

i 

deg 

Rs 

kpc 

PO 

stars pixeU^ kpc“^ 

PBG 

All 

- 

129.51 ±0.17(±1.08) 

-38.14 ±0.08(±1.59) 

1.09 ± 0.01(±0.02) 

26.41 ±0.14(±0.59) 

0.80±0.03(±0.67) 

Young 

0 - 4 

125.93 ±0.20^0.09) 

-44.19 ±0.14^1-80) 

0.72±O.OiLo.OO) 

22.86 ±0.28(±1-33) 

9.06±0.18(±0.00) 

Old 

4 - 13 

127.40 1 tl. 02 L 0 . 59 ) 

-32.94 ±0.39Li-25) 

1.41 ± O.OlLo.OO) 

0.89 ± 0.02(±0.05) 

2.08 ± 0.04^0.01) 


Table 3. A summary of the best-fit disk models parameters for each stellar population. The uncertainties are the 3cr confidence level 
that arises from the MCMC analysis. The uncertainties in parenthesis are obtained from the difference between the fit with and without 
correcting the magnitudes by the heliocentric distance to each point of the LMC disk. NOTE: our definition of inclination (i) has the 
opposite sign than what is typically found in literature. 


to the Galactocentric distance of the LMC (^lmc)- 

Myrw{d < o^LMc) = 24.5l6;4 X Mlmc (5) 

where the Sun is assumed to be at a distance of d© = 8.0 
kpc from the MW centre. This yields a o^lmc of 49.4 zb 2.1 
kpc. The uncertainty in this value was considered for the 
result in Equation 5. The distance to the Sun is taken as a 
compromise between two recent determinations from Eisen- 
hauer et al. (2005) and Gillessen et al. (2009). The adoption 
of d© = 8.5 kpc would increase o^lmc by ~ 0.03%, leading to 
an insignificant increase in Rt. Thus, we choose to disregard 
the uncertainty in d©. 

To further support our claim that there are very few 
LMC stars beyond 13 kpc we show in Eigure 10 the Hess di¬ 
agrams at different bins of angular separation from the LMC 
centre. The upper and lower bound of each bin is indicated 
in each panel. We apply a simple decontamination algorithm 
in order to remove the contribution of MW stars. The de¬ 
contamination was done by selecting a region with angular 
distance greater than 0 — 20°. A Hess diagram of this region 
was constructed. These Hess diagrams were area weighted 
and subtracted from the Hess diagram of each angular sep¬ 
aration bin shown in Eigure 10. This process assumes that 
the contribution of LMC stars at angular distances larger 
than 20° is negligible. We again assume that the MW stel¬ 
lar population varies very little throughout the SPT-E. 

In Eigure 10 we note, from visual inspection, that the 
LMC population in the bottom central panel is indistin¬ 
guishable from the noise in the Hess diagram. This panel 
corresponds to angular separations between 15° and 16°. 
This range corresponds to distances from 10.5 and 11.2 kpc 
along the disk plane, which are intermediate between the 
truncations radius of the old and young population. This 
does not mean the complete absence of LMC stars. However, 
it shows that at this range of distances the LMC populations 
is too low in number to be distinguishable by eye from the 
MW foreground stars. 

With the aid of the GALPY^ (Bovy 2010) suite for galac¬ 
tic dynamics we compute a three component MW poten¬ 
tial with a NEW halo (Navarro, Erenk & White 1996), 
a Miyamoto-Nagai disk (Miyamoto Sz Nagai 1975), and a 
Hernquist bulge (Hernquist 1990) following the same recipe 
as in Bovy et al. (2012). This potential was normalized such 
as to yield a Solar circular velocity of 220 km/s at the Galac¬ 
tocentric distance of 8 kpc. We obtain that the enclosed MW 
mass inside a sphere centred in the MW with radius dLMC 


^ http://github.com/jobovy/galpy 


is 5.6 X lO^^M©. Using the result from Equation 5 and prop¬ 
agating the uncertainties in Rt we obtain an LMC mass of 
Mlmc = 2.3l“;« x 

According to the set of simulations by Kallivayalil et al. 
(2013), the combination of MW and LMC masses found in 
this work would favor a scenario where the Magellanic Sys¬ 
tem is in its first pericentric passage. However, it is not clear 
if there is enough time for the LMC to develop a truncation 
radius without finishing a complete orbit around the MW 
and if this truncation would affect its luminous component. 
We would like to stress that the calculations presented in 
our work are subject to many uncertainties, especially with 
respect to our knowledge of the MW mass profile at large 
distances. 


5 THE LMC RED CLUMP AS A DISTANCE 
ESTIMATOR 

The peak magnitude of the RC is used widely in astronomy 
as a standard candle to determine distances of old and in¬ 
termediate age star clusters. RC stars are He-burning stars 
of mass M < ‘2Mq', they develop an electron degenerate 
core after the main sequence, and as a consequence have to 
increase their core up to a critical mass of ~ 0.46M© be¬ 
fore Helium can be ignited. The almost-constancy of their 
Helium-core masses determines that they all share similar 
(but not constant) luminosities in their central helium burn¬ 
ing phase. 

Using the RC peak magnitude as a standard candle we 
define the distance modulus /tq with the following equation: 

Ato = -Ax- AMf ^ (6) 

where Ax is the extinction, is the absolute magnitude 

of the RC peak, and AMx^ is a correction in the absolute 
magnitude due to population mixing effects. Here A repre¬ 
sents the observational filters (A = g, r, z, z). 

If one is to determine /xo, the terms Ax and AMx^ 
must be inferred. Here we choose to correct the magnitudes 
for extinction using the Schlegel, Einkbeiner & Davis (1998) 
dust maps. Eor the time being we will assume that AMx^ 
is zero. This assumption will be addressed in more detail 
later. 

To probe /tq as a function of position in the LMC, we 
subdivided the sky in HEALPIX pixels using the pixelization 
scheme with Nside = 128. This gives a pixel area of ~ 0.21 
deg^. We choose a larger pixel than the previous section since 
the RC is much less populated than the Main Sequence. To 
measure the peak apparent magnitude of the RC in a given 
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Figure 6. Left panel: map showing the number of stars from the all sample in grayscale. The stars used to build this map are those that 
fall inside a simple color cut of —0.5 < g — r < 1.0 and g < 23. The solid isodensity contours show the best exponential disk model. The 
contours start at log]^Q(A/^) = 3.3 and progress in steps of 0.3 dex. Right panels: the marginalization for the different disk model pairs of 
parameters. The solid blue crosshair shows the position of the best solution. The contours show the 1, 2, and 3cr confidence levels. The 
histograms in the diagonal panels show the marginalization over each single parameter. Again the blue line shows the point of maximum 
likelihood. 



Figure 7. Similar to Figure 6 but for the young stellar population. The isodensity contours start at log]^Q(A/’) = 2.9 and progress in 
steps of 0.4 dex. 
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Figure 8. Similar to Figure 6 but for the old stellar population. The isodensity contours start at log]^Q(A/’) = 2.5 and progress in steps 
of 0.35 dex. In the left panel we also show the truncation radius (Rt) position as the dashed contour, and the dotted contours show the 
uncertainty in Rt. 


passband (rrix^) we compute the number count of stars as a 
function of magnitude {N{mx)) for stars that have colours 
and magnitudes limited by the dashed boxes in Figure 4. The 
CMD region occupied by RC stars was selected visually. 

To build N(mx) the bin size is chosen according to the 
method described in Knuth (2006), which is based on op¬ 
timization of a Bayesian fitness function across fixed-width 
bins. This avoids issues that may arise from oversampling or 
undersampling the number of bins. The values of N{mx) are 
fitted by means of a non-linear least square algorithm using 
a second order polynomial plus a Gaussian. This function is 
given by the following equation: 


^(^a) = a + bmx + cm\ + dexp 




(7) 


where a, 5, c are the coefficients of the polynomial, d is the 
normalization of the Gaussian and ax is the standard devi¬ 
ation of the magnitudes around the RG peak. 

The uncertainty on rux^ is taken from the covariance 
matrix of the least square fit. 

It is convenient to define the heliocentric distance to 
points in the LMG disk as a function of the so called line 
of maximum distance gradient. This line connects points 
on the disk plane with the most rapidly varying distance 
from us, hence its name. The distance along this line is the 
deprojected distance between a point in the disk to the line 
of nodes. The distances along this line are given by the y 
component of Equation A3 of Weinberg & Nikolaev (2001). 
The line of maximum distance gradient has a PA given by 
0 (Equation 2) and it is oriented approximately in the NE- 
SW direction. Showing the heliocentric distance to points 
in the LMG disk as a function of the position in the line of 


maximum distance gradient is the equivalent of showing an 
edge-on view of the LMG. 

To infer the absolute magnitude of the RG (M^^) we 
could adopt the prediction from synthetic stellar populations 
based on stellar evolution models. However, simulations us¬ 
ing PARS EG (Bressan et al. 2012) models agree well with 
the disk models found in this work. A small magnitude offset 
can be seen on Eigure 11, where the solid black line shows 
the best-ht disk model and the gray points are the inferred 
distance modulus based on theoretical values of Mx^. 

Despite the small offset in the theoretical value of Mx^ 
we chose to determine Mx^ by matching the expected dis¬ 
tance modulus from our best-ht disk model to the observed 
values of the RG peak magnitude. This matching was done 
using stars that are between 3 and 4 kpc along the LMG 
line of maximum distance gradient. We compute the me¬ 
dian value of mx in this distance range and subtract from 
that the heliocentric distance expected from the disk model 
at 3.5 kpc from the LMG centre. This offset is determined for 
ah passbands. The value obtained is then used to compute 
the RG-based distances consistently with our disk model. 
These distances are shown as the black points in Eigure 11 
and the error bars are propagated from the RG peak ht. 

In Eigure 12 we show the distributions of rrix^ on the 
sky, for A = g, r, i, z. The lines show the direction of maxi¬ 
mum gradient and the line of nodes obtained from our best 
disk model ht using the all stellar population. We notice 
that the entire sampled region in this work is on the same 
side (near side) of the LMG disk major axis. There is a 
global trend in the sense that for ah passbands the nearest 
points to us are located in the North-East edge of the LMG. 
Another remarkable feature present in ah maps is that the 
North edge of the LMG is systematically more distant than 


© 0000 RAS, MNRAS 000 , 1-18 


























































12 E. Balbinot et al. 



0.5 0.6 0.7 0.8 0.9 1.0 1.1 1.2 

logio(R/kpc) 


Figure 9. The black circles show the averaged number of stars 
from the all sample per HEALpix pixel in different bins of dis¬ 
tances along the LMC disk. The triangles (crosses) show the same 
but for the old {young) sample. The dashed, dot-dashed, and dot¬ 
ted lines show the best-fit disk model for the all, old, and young 
populations, including the term pBGf{(^i ^) that accounts for the 
contamination of MW stars. In the inset plot we show the outer 
tail of the old profile. The solid line shows the MW contribution 
and the shaded region shows its uncertainty, propagated from 
the uncertainty in pbg- The errobar represents the Poission un¬ 
certainty on the number of counts per bin. 



Qi-Qi* Qi-V Qi* Qi* Qi* V 

g-r g-r g-r 


Figure 10. Hess diagrams in — r vs. g for different angular 
distance ranges (indicated in the top left of each panel). The 
contribution from MW stars has been removed as described in 
the text. 


what is expected for the disk model. This effect is noticeable 
in all passbands. 

The discrepancies observed in the maps described above 
may be due to population mixing in RC stars. The RC lu¬ 
minosity depends slightly on the age and metahicity of the 
stars (Girardi 1999; Girardi & Salaris 2001). This fact makes 
the RG of a stellar population older than ~ 3 Gyr dimmer 
as a function of age. The exact amount of dimming depends 
on the metahicity and on the passband used. A detailed 
discussion about how the RG magnitude of simple stellar 
populations (SSP) varies with age and metahicity may be 
found in Girardi & Salaris (2001). 

In the case of the LMG we would like to know how the 
RG magnitude changes as a function of the SFH, therefore 
mapping at different positions in the LMG. However, 

it is known that the LMG SFH is very complex and varies 
spatially (Meschin et al. 2014), therefore rendering this a 
very difficult task. 

To assess the amplitude of the populations mixing effect 
over the RG magnitudes, we conducted a series of synthetic 
stellar population experiments. These simulations assume 
that the LMG follows an age-metahicity relation given by 
Piatti & Geisler (2013) with a spread of 0.15 dex in [Fe/H]. 
We also assume that stars in the LMG follow a Kroupa Ini¬ 
tial Mass Function (IMF) (Kroupa 2001). To generate syn¬ 
thetic magnitudes for a set of simulated stars we adopt the 
PARSEG (Bressan et al. 2012) isochrones in the DES pho¬ 
tometric system. Our grid of models has a step of 0.01 in 
log^Q(a^e/yr) and 0.0002 in Z in the range [0.0002,0.001] 
and 0.001 for Z in the range [0.001,0.020]. Photometric un¬ 
certainties were included according to Equation 1 and Table 
2. The simulations were made in the framework of the open- 
source code genGMD^ 

We adopt three SEHs that are modelled according to 
what has been found by Meschin et al. (2014) for the outer 
regions of the LMG (red squares in Figure 12). The simu¬ 
lated SEHs contain two events of star formation that were 
modelled as Gaussian peaks in SFR. The young peak is cen¬ 
tered on an age of 2.5 Gyr with a width of 1 Gyr while the 
second is centered on 9 Gyr with a width of 1.5 Gyr. The 
hrst adopted SFH gives equal amplitude to each star form¬ 
ing event, the second gives twice the amplitude to the older 
event, and the third gives twice the amplitude to the younger 
event. The largest difference for the RG peak magnitude was 
observed between the models with asymmetric peaks in the 
SFH. We use this difference to estimate the maximum value 
for AMx^ which was found to be 0.06, 0.06, 0.07, and 0.11 
for g, r, z, z respectively. The dashed lines in Figure 11 show 
the maximum expected deviations due to populations mix¬ 
ing effects. 

We notice that the RG distance variations are very simi¬ 
lar to the one expected from the best-ht disk model. However 
there are very discrepant points that the disk model cannot 
account for, even when population mixing effects are con¬ 
sidered. We also notice that the distance moduli determined 
using different passbands are consistent. 

From the RG distributions we may also extract informa¬ 
tion about the thickness of the LMG disk. This information 
is embedded in ax shown in Equation 7. However, this quan- 


^ https://github.com/balbinot/genCMD 


© 0000 RAS, MNRAS 000 , 1-18 







The LMC outer stellar populations and geometry 13 






Figure 11. The distance modulus as a function of the distance along the maximum gradient line NE-SW). The black circles are the 
measured distance modulus in each HEALpix pixel, as describe in the text. The error bars shows the 1 CTm for each measured point. 
The solid line shows the behaviour expected from the best-fit disk model for the all population. The dashed line shows the maximum 
magnitude variation expected from populations mixing effects. The gray points are the distance modulus inferred using the theoretical 
RC absolute magnitude. Panels a,b,c, and d show the distance modulus determination using the filters g, r, i, and z respectively. 


tity is also affected by the intrinsic scatter of the RC and 
by photometric uncertainty. The RC simulations described 
above incorporate this intrinsic scatter (ai) convolved with 
the photometric one (aphot)- In our simulations we measure 
this quantity (and we find that this quantity 
has a mean value of 0.11, 0.12, 0.13, 0.13 for ^,r, z, and z 
respectively. Finally the last contributions comes from the 
depth of the disk along the line of sight itself {cf depth)- Thus, 
we may write: 


2 _ 2 2 2 
^depth ^phot - 


(8) 


In Figure 13 we show a depth as a function of the dis¬ 


tance along the LMC disk. The error bars are taken from 
the covariance matrix obtained on the least-square fit of the 
RC distribution as given by Equation 7. We took averages 
of crdepth in bins of 0.5 kpc, which are shown as the black 
circles. The error bar is the standard deviation of the values 
of cr depth in each distance bin. For distances less than ~ 4.5 
kpc we observe a constant depth of ~ 0.08 mag (1.8 kpc) 
measured in r, z, and z while the g band yields a value of 
0.12 mag (2.8 kpc). The depth in r, z, and z passbands tends 
to increase by ~ 0.02 mag (~ 0.5 kpc) towards the edge of 
the LMC. The slight difference in the g band depth might 
be due to the underestimate of the intrinsic RC magnitude 
spread. This effect may be due to changes in the metallicity 
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Figure 12. Gnomonic projection of Ngide = 128 maps showing the distance modulus measured in each the and z passband. On 

every panel the solid line shows the direction of the maximum distance gradient (6) and the dashed line the direction of the line of nodes 
(^o) expected from the best-fit disk model for the all population. The red squares are the fields studied by Meschin et al. (2014). The 
red circle is the LMC centre adopted in this work. Panels a, b, c, and d show the distance modulus determination using the filters g, r, 
z, and z respectively. 


of the LMC field populations torwards its edge (Majewski 
et al. 2009), leading to the formation of a Horizontal Branch 
(HB) instead of a RC. The g band is more affected since it 
is more sensitive to hotter stars. 


6 DISCUSSION AND SUMMARY 

To the best of our knowledge, this paper is the first to show a 
clear distinction in the structure of the LMC disk for stellar 
populations of different ages. Meschin et al. (2014) report 
signs of an age gradient, but did not derive the structural 
properties for different age components of the LMC. Previ¬ 
ous studies, based on large spectroscopic samples, reported 
a clear kinematic distinction between different stellar types 
and attributed that to an age difference in their samples (van 
der Marel, Kallivayalil & Besla 2009, and references therein). 
Here we confirm this distinction based only on stellar pho¬ 
tometry alone. 

The disk models obtained in this work are in agreement 
with what has been reported in the literature (e.g. Rubele 
et al. 2012; Saha et al. 2010; Nikolaev et al. 2004; Weinberg 


& Nikolaev 2001). We observe a significant difference be¬ 
tween disk models fitted using stars younger and older than 
4 Gyr. The most striking difference is found in the scale 
length (Rs), which suggests a star formation that follows 
the outside-in paradigm. The summary of the disk models 
parameters can be found in Table 3. 

The young population is poorly described by an ex¬ 
ponential disk model. It presents a prominent feature that 
strongly deviates from an exponential profile truncated at 

kpc. The young and old LMC disk have a clear dis¬ 
tinction in the level of substructures. The younger stars are 
known to have peculiar spiral arms that have most likely 
been formed in the last SMC-LMC encounter (Staveley- 
Smith et al. 2003; Olsen & Massey 2007; Bekki 2009). This 
could explain why the young disk cannot be well described 
by our disk model. 

On the other hand, the old disk extends over ~ 13 scale 
lengths. Our results for the scale length and truncation ra¬ 
dius are similar to what Saha et al. (2010) found. The au¬ 
thors report Rs ^ 1.15 and a truncation radius that is ~12 
times this value, which is very similar to our determination. 
If we assume that the truncation radius of the old disk is 
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Figure 13. Disk depth as a function of distance to the LMC centre. The gray error bars show the Cdepth measured for each HEALpix 
pixel shown on Figure 12. The dark circles are the average over 0.5 kpc and the error bar the standard deviation. Panels a,b,c, and d 
are for determinations using the filters and 2; respectively. 


caused by the tidal field of the MW, we obtain a LMC mass 
that is Mlmc — 2.3lQ g x This mass is strongly 

dependent on the efficiency with which the MW potential 
tidally truncates the LMC stellar distribution. Even if the 
Clouds are in their first pericentric passage, their orbit is 
predicted to be past the pericenter, which could allow the 
formation of a tidal limit (Besla, private communication). 
Also, (Gan et al. 2010) showed that subhalos in their first 
pericentric passage may undergo substantial mass loss. We 
also point out that in the case where the MW tidal field 
efficiently strips the dark matter and luminous components, 
our estimate of the LMC mass is reliable. However, in the 
case where dark matter is still present in an extended halo, 
our estimate of the LMC mass sets only a lower limit. The 
absolute mass value is also dependent on the MW mass pro¬ 


file, but the relative LMC to MW encircled mass is not. New 
more precise determinations of the MW mass in its outer re¬ 
gions will help to more accurately constrain the LMC mass. 

The dynamical mass determination by van der Marel & 
Kallivayalil (2014) suggests Mlmc{R < 8.7 kpc) = (1.7 zb 
0.7) X 10^^ Mq that can be extrapolated to the mass inside 
our truncation radius of Mlmc{R < Rt) = (3.5 zb 1.4) x 
10^°M 0 assuming a flat rotation curve. This value is within 
the error bars of our mass determination, suggesting that 
the truncation radius found in this work has a tidal origin. 

Through the fit of the RC peak magnitude we observe 
that a few regions coherently grouped in the North of the 
LMC are systematically more distant than what is expected 
by an inclined circular disk model. While Olsen & Salyk 
(2002) found that regions of the South-East LMC are sys- 
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tematically closer, we find that the Northern LMC disk be¬ 
haves in the opposite sense. This is the classic case for a 
galaxy with warped disk. We demonstrate, through the use 
of synthetic stellar population models, that the offsets from 
the disk model are larger than what is expected from popu¬ 
lation mixing effects. 

Using the same synthetic stellar populations quoted 
above we are able to compute the theoretical magnitude 
scatter of RC stars. This allows us to disentangle the con¬ 
tribution of the intrinsic scatter from the one caused by the 
disk thickness. We measure a thickness ranging from 1.8 to 
2.8 kpc in the inner parts of the LMC and which increases 
by 0.5 kpc in its outer parts. If we assume that the disk 
follows an exponential profile in height we obtain that the 
scale height ranges from 1.3 to 1.9 kpc, increasing by 0.3 kpc 
in the outer parts. This slight increase towards the outskirts 
can be interpreted as the flaring of the disk, van der Marel 
& Cioni (2001) argue that this kind of behaviour is expected 
if the LMC disk is tidally perturbed by the MW potential. 
The disk thickening may also be explained solely by LMC 
dynamical models (Besla et al. 2012), where the old stellar 
component is expected to form a dynamically hotter com¬ 
ponent. The fact that we only observe it on the edge of the 
LMC may be due to the steep truncation of the young disk, 
allowing us to measure a thickness that is dominated by the 
old population only in the LMC outskirts. 

The analysis presented in this work points to a scenario 
where there are two distinct disk components in the LMC. 
One is composed of old stars (> 4 Gyr) and possesses a 
smooth, extended profile out to dozens of scale lengths. The 
second component is composed of younger stars (< 4 Gyr) 
and appears to be pertubed and much less extended. We find 
signs of warp and flare towards the outskirts of the LMC, 
both of which have been reported previously (Olsen & Salyk 
2002; Subramaniam & Subramanian 2009). 

An alternative scenario for the disk warp and flare is the 
existence of a hot halo composed of old stars, while the disk 
could contain both populations. Because the old population 
is more spatially extended, this scenario would explain the 
outward increase in the thickness of the RC as the young disk 
fades, therefore accounting for the “flare” signal. Besides 
the flaring^ the old spheroidal component would lead to a 
systematic variation on the heliocentric distances across the 
LMC in the sense of bringing them closer to the heliocentric 
distance of the LMC centre on both sides of the galaxy, also 
mimicking a warp effect. 

The fact that the old disk is well-fit by a single expo¬ 
nential disk favours the scenario where there is no detectable 
stellar halo and older stars form a thicker disk. Also the fact 
that the hot halo has only been possibly detected in the cen¬ 
tral regions of the LMC (Minniti et al. 2003) raises doubts 
about the existence of such a stellar component. To inves¬ 
tigate the existence of a halo component, a more detailed 
study must be carried out in the LMC outskirts in order to 
isolate its extremely sparse stellar population. This kind of 
study should employ more sophisticated CMD decontami¬ 
nation methods such as the ones used to detect streams and 
tidal tails (Rockosi et al. 2002). 

The upcoming years promise to greatly increase our un¬ 
derstanding of the Magellanic system. With the data already 
observed during the DES SV campaign, we have the ability 
to study a vast number of LMC star clusters and their 3- 


dimensional distribution in the galaxy. The field LMC popu¬ 
lation may also be used to infer a very detailed, and spatially 
dependent SFH. 

Upcoming DES observations will cover a large portion 
of the Magellanic Stream allowing for the discovery and 
characterization of its stellar component, if it exists. Addi¬ 
tionally, DES will continue to map the outskirts of the LMC 
and SMC. This novel dataset will certainly reveal some of 
the many puzzles of the Magellanic system. 
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